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Abstract 



We investigate spontaneous interlayer phase coherence and the occurrence 
of the quantum Hall effect in triple-layer electron systems. Our work is based 
on a simple tight-binding model that greatly facilitates calculations and whose 
accuracy is verified by comparison with recent experiments. By calculat- 
ing the ground state in an unrestricted Hartree-Fock approximation and the 
collective-mode spectrum in a time-dependent Hartree-Fock approximation, 
we construct a phase diagram delimiting regions in the parameter space of 
the model where the integer quantum Hall effect occurs in the absence of 
interlayer tunneling. 

PACS numbers: 73.20.Dx, 75.10.-b, 64.60.Cn 
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I. INTRODUCTION 



Interactions among particles play an especially important role in two-dimensional (2D) 
systems in the quantum Hall (QH) regime because the kinetic energy is quenched. Kinetic 
energy eigenstates bunch into quantized Landau levels with macroscopic degeneracy N$ = 
B±A/§q. (Here A is the cross-sectional area of the system, B± is the magnetic field strength, 
and $o — hc/e is the magnetic flux quantum.) The fractional QH effect (QHE),0 which 
occurs when an orbitally degenerate Landau level has a fractional filling factor z/p = Ng/N^, 
where N e is the total number of electrons, arises from strongly correlated states produced 
entirely by interactions]!! 

Interactions can play an important role even at integer values of z/r, if at low energies the 
system has additional degrees of freedom. An important example occurs for an isolated 2D 
electron layer at z/r = 1, when the Zeeman energy is so small that the spins are not completely 
frozen. It turns out in this case that the ground state is completely spin polarized, and that 
the energy gap for charged excitations (which gives the QHE) is finite, even in the limit 
of vanishing Zeeman splitting. However, there is strong experimental evidence that for z/p 
close, but not equal, to 1, the ground state contains a large number of flipped spins.@ This 
property of the single-layer system, which was anticipated theoretically,i~@ is best understood 
by recognizing that the ground state at z/r = 1 is a 2D ferromagnet with spontaneous spin- 
polarization. 2D ferromagnets have stable finite-energy topologically charged spin-texture 
excitations, commonly known as skyrmions.i A unique aspect of these QH ferromagnets, 
first appreciated by Sondhi et al.p is that their skyrmions also carry a unit electrical charge. 
It is the presence of these unusual charged objects with many reversed spins in their interior 
that is responsible for the rapid decline in the spin polarization that occurs as z/r moves 
away from 1. Many aspects of the physics are very similariHl! when the additional degree of 
freedom comes from a second electron layer, rather than from the two spin states available 
to spin- 1/2 particles. The role played by the Zeeman energy is taken over in this case by 
the interlayer hopping amplitude t, the broken symmetry is spontaneous interlayer phase 
coherence rather than spin magnetization, and a QHE can occur0lll at z/r = 1, even when 
t = 0. The combination of the integer QHE and spontaneous broken symmetry in the ground 
state gives rise to charged order-parameter textures and other new physics.EI 

In this paper, we consider the case of triple- layer electron systems (TLESs). Our work is 
motivated primarily by recent progress in fabricating high-mobility electron systems, which 
has made it possible to study these systems experimentallyt^rEJ Triple-layer systems in 
strong magnetic fields have been studied previously,!]! but the possibility of QHEs associ- 
ated with spontaneously broken symmetries,0 on which we focus here, was not explicitly 
addressed. In Sec. II of this paper, we introduce the tight-binding model that we use to 
describe TLESs in terms of a small number of parameters. We make some estimates of 
the size of the model parameters in Sec. Ill, and discuss some electrostatic considerations 
that are very important in interpreting experiments. We test our model and demonstrate 
the possibility of determining its parameters by comparing with some recent experiments 
at weak magnetic fields. In Sec. IV, we discuss unrestricted Hartree-Fock approximations 
(HFAs) for the ground state of TLESs in a strong magnetic field. The HF wave functions 
allow for the possibility of spontaneous interlayer phase coherence and are generalizations 
of those proposed previously^ for double- layer systems at z/ T = 1: 
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l*> = n(^x + ^ 2 cL)|0>. (1.1) 
x 

Here X is the guiding-center label for orbital states within a Landau level in the Landau 
gauge, which we use throughout this paper. This many-particle wave function exhibits 
interlayer coherence, because all electrons occupy states that are a linear superposition of 
layer j — 1 and layer j = 2. When <pi = <ft 2 , Eq. ( |1 . 1| ) represents a full Landau level formed 
from the symmetric combination of isolated layer states, and is evidently the exact ground 
state for a system of noninteracting electrons, when tunneling between the layers is described 
in a tight-binding model. It turns out that it is also the exact ground state in the presence 
of repulsive interactions, even without tunneling, when the interactions are independent of 
the layer index j (i.e., in the limit that the layer separation d — > 0). When the interlayer 
Coulomb interactions are different from the intralayer interactions (d > 0), Eq. (|1 . 1| ) is 
still a good variational wave function because phase coherence guarantees good interlayer 
electronic correlations and thereby lowers the interlayer Coulomb interaction energyi Only 
when the layers are widely separated does Eq. ( |l.lj ) become a poor variational wave function; 
eventually, it becomes much more important to have good intralayer correlations than to 
have good interlayer correlations. In this case, the ground state no longer has spontaneous 
interlayer correlations and the QHE will not generally occur at z/p = 1 in the absence of 
interlayer tunneling. (An exception occurs when the electrostatic environment is consistent 
with having QHEs in which interlayer correlations play no role, for example by having 
v — 1/3 in each layer.) In Sec. V, we evaluate the collective-mode dispersion of a TLES in 
the time-dependent HF approximation (TDHFA). In the absence of interlayer tunneling, the 
energy in the double-layer case is independent of the global phases <pj. The U(l) symmetry 
associated with <p 12 = (pi — (p 2 is broken in the ground state of the double-layer system, 
giving rise to a gap less Goldstone mode associated with slow spatial variations of 0^2 -EM 
The energy cost of the spatial variations of cf>i 2 (X) is due to the loss of interlayer Coulomb 
exchange energy. We find that in triple-layer systems, two Goldstone modes occur, associated 
with the two independent relative phases. We use the stability of the collective modes as 
an indication that the HF variational wave function is still good, and use this criterion to 
map out the region in the model's parameter space where we expect a triple-layer integer 
QHE to occur with spontaneous interlayer phase coherence in the ground state. Section VI 
contains some brief concluding remarks. A future work will discuss charged order-parameter 
textures and the effect of tilted magnetic fields in triple-layer systems, using a field-theoretic 
approach. 



II. MODEL FOR TRIPLE-LAYER SYSTEMS 



Triple-layer QH systems have been realized experimentally by Shayegan and co- 
workers.tjUa In order to observe the QHE, the mobility of these samples needs to be very 
high, which generally requires that they be remotely doped. The simplest possible theory of 
the triple-layer system is one that regards it as a macroscopic metal, and therefore requires 
that the volume between the right and left layers be an equipotential. In such a theory, 
the Poisson equation allows no charge in the central layer; repulsive interactions between 
the electrons cause them to migrate to the left and right layers. Obviously, this theory is 
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incomplete, but its indication that electrons tend to avoid the middle layer is telling, and 
this tendency must be countered if true triple-layer systems are to be realized. In the TLESs 
grown by Shayegan's group, the middle layer is wider than the outside layers and therefore 
has a smaller size-quantization energy. 

The Shayegan group has demonstrated@~El that at zero magnetic field, the partitioning 
of electron density between the three layers as a function of gate voltage can be accurately 
rendered using a three-dimensional (3D) density-functional based independent-electron ap- 
proximation. In order to describe the many-body physics of these systems, which is essential 
at strong magnetic fields, we require a relatively simple model for the growth-direction ( in- 
direction) spatial degree of freedom. We proceed by generalizing the approach commonly 
used for double-layer 2D electron systems. We assume that only the lowest electric sub- 
band is important in each 2D layer, and use a tight-binding description for the ^-direction 
degree of freedom, with tunneling amplitude t between neighboring layers, and no direct 
tunneling between left and right layers. (Remote tunneling can be readily incorporated 
in the model, if the experimental situation warrants introducing this complication.) For 
the calculations reported here, we neglect the finite width of the subband wave functions 
in each layer, but these can easily be modeled if necessary in specific systems by adding 
form-factor corrections^ to the effective electron-electron interaction. We assume that the 
subband energy in the middle layer (relative to the local electrostatic potential) differs from 
the subband energies in the side layers (also relative to the local electrostatic potential) by 

In experiments on TLESs, it is extremely useful to be able to manipulate charges in the 
layers with both a front (F) gate, which we take to be closest to the left layer, and a back 
(B) gate, which we take to be closest to the right layer. We parameterize these gate voltages 
in terms of neutralizing charge densities p a (a =F,B), by defining 

V aG = —p a D aG -V$!, (2.1) 
e o 

where e ~ 13 is the dielectric constant for GaAs and D aG is the distance from the closest 
layer in the system to gate a. (For the systems fabricated by the Shayegan group, -Dfg ~ 
0.45 /im and -Dbg ~ 0.45 mm.E3) The last term represents an offset voltage. This model for 
the triple-layer system is illustrated schematically in Fig. [I]. 

To test the appropriateness of such a model, and to determine model parameters for 
the experimental system with which we will later compare our strong magnetic-field theory, 
we have calculated the dependence of the state of the triple-layer system on the front-gate 
voltage (Vpc). For this purpose, we have used a 2D version of the 3D local-density-functional 
(LDF) theory.il The Kohn-Sham single-particle equations separate and yield three two- 
dimensional free-electron bands with minima at subband energies that are determined by 
solving the three-site discrete Schrodinger equation for the z-direction degree of freedom. In 
this equation, neighboring layers are coupled by the tunneling matrix element —t and the 
LDF site energies are given (up to a common constant term) by 

2ne 2 d 

£i = [(ni - n 3 ) - (p F -pb)] + fhcM, 

eo 

2iie 2 d , . 

e2 = e& H n 2 + // xc (n 2 J, 

eo 
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£3 = [(n 3 - Tii) - (p B - Pp)] + H xc (n 3 ). (2.2) 

The "bare" middle-well on-site energy e& represents the difference in the size-quantization 
energy of the middle layer, relative to that of the side layers. It is negative (attractive) 
when the middle well is wider than the side wells. The Hartree terms are proportional to 
27ie 2 d/e . The areal densities of the left, middle, and right layers are denoted by n\, n 2 , 
and 7^3, respectively. The exchange contribution to the exchange-correlation potential /i xc 
of local-density-functional theory is fi x (n) = d[ne x (n)]/dn = — e 2 y / 8n/7r/eo, where e x (n) is 
the contribution to the exchange energy per particle in a 2D electron gas of uniform density 
n@ We have not included correlation-energy contributions to /i xc , but this could be added if 
desired. The subband energies E\ are obtained by diagonalizing the 3x3 LDF Hamiltonian 
matrix. The density in layer j is given by 

n s = E# A |af | a , ( 2 - 3 ) 

A=l 

where 

N x = (E F - E x )u Q(E F - E x ) (2.4) 

is the areal-density contribution from the Ath subband, is the amplitude of the Ath 
subband wave function in layer j, and Vq = m* /irh 2 is the 2D electron-gas density of states 
(m* ~ 0.07m e is the band effective mass for GaAs). Equations (|2.2|) and Q2.3|) have to 
be solved self-consistently, with the Fermi energy chosen so that J2\N\ = Nt = Pf + Pb, 
consistent with overall charge neutrality. 



III. MODEL PARAMETERS 

Rough estimates of the size of the model parameters may be obtained from simple argu- 
ments. In the limit of infinitely strong barriers separating the layers, the size-quantization 
energies are Tr 2 h 2 /2m*w 2 , where Wj is the width of the jth quantum well. For side- well widths 
of 15.4 nm and a middle-well width of 18.8 nm, the difference in the size-quantization en- 
ergies gives ~ —8 meV for the TLES of Shayegan and co-workers that we study here. 
We note that the next group of three electric subbands is higher in energy by roughly 
w 2 ~ 50 meV, and may be neglected in studying the ground state and low- 
energy excitations of the TLES. The tunneling energy may be crudely estimated from a 
semiclassical (WKB) argument, by equating 2t/H with the rate at which the wave-function 
amplitude of an electron leaks out of its well by tunneling through the confining barrier: 

2 4 ~ -H-e-"* = \*e-* (3.1) 
h 2w h vr K ' 



where v/2w is the frequency with which an electron of average velocity v = ^2E /m* in 
a well of width w ~ 17 nm hits the side of the well, Eq ~ ti 2 Ti 2 /2m*w 2 ~ 20 meV is the 
energy of the confined electron, hn = y / 2m*(Vo — Eq) ~ y/2m*Vo, whereil Vq ~ 1 eV is the 
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barrier height, and b ~ 1.3 nm is the width of the barrier. From this estimate, the tunneling 
energy t ~ 0.5 meV, and is thus expected to be an order of magnitude smaller than the 
on-site energy e&. We numerically solved the one-dimensional Schrodinger equation for the 
triple quantum-well potential shown in Fig. [I] in an effective-mass approximation to obtain 
the bound-state energies. The three lowest energies determine e b and t through Eqs. ( p.5| ) 
and ( |3.6|) , and give e b ~ —4.7 meV and t w 0.5 meV. We also find an energy separation of 
34 meV to the next group of three electric subbands. 

Electrostatic energies tend to be larger than planar kinetic energies, tunneling energies, 
and exchange-correlation energies. For the densities shown in Fig. electrostatic energies 
are typically about four times larger than planar kinetic energies, about six times larger than 
the exchange energy, and can be on the order of 100 times larger than the tunneling energies, 
depending on the barrier widths. As a result, Fig. ^| can be understood qualitatively in an 
electrostatic approximation, where the electronic charge resides entirely in the layer with 
the lowest (Hartree plus "bare" on-site) site energy, unless two or more layers are placed in 
equilibrium by having the same energies. For 4ire 2 dpF/eo < e& < 0, it follows that in this 
approximation, all electrons will occupy the right layer, and n 3 = Pb+Pf- With increasing 
Pf, charge will be added first to the middle layer, until e b + 4ne 2 dpF/eo « 0. As pf is 
increased further, all charge is added to the left layer. These considerations provide two 
straightforward measurements of e b . First, at the front-gate voltage V3 ~ —0.3 V where the 
third subband first becomes occupied, electrostatic considerations show that 

e b « -^^N 2 (V 3 ) « -4.7 meV (3.2) 

for ^2(^3) ~ 3.5 x 10 10 cm -2 and d = 18.4 nm. Second, for fixed back-gate voltage, the 
difference between V3 and the value of the front-gate voltage (V 2 ~ —0.6 V) at which the 
second subband first becomes occupied also measures e b : 

e b « --7^(^ 3 - V 2 ) « (y 3 - V 2 )-/- ~ -6 meV (3.3) 

Z L>FG e « ^FG 

for V3 — V 2 ~ 0.3 V. In Eq. ( [3. 3D , we have used the fact that 1/Dfg is proportional to the 
front-gate capacitance per unit area, d(— eN T ) / dVFG- 

When the dependence of the chemical potential on the density in each layer is taken into 
account, there is a small correction to the electrostatic result. This correction is measurable, 
and has been exploited by Eisenstein and co-workersSi! in the double-layer case to measure 
the compressibility of the electron-gas systems within each layer. For example,^ features 
occur in the charge- density distribution when the density in one of the layers is very small, 
which reflect the diverging (negative) compressibility in the low-density limit of an electron 
gas. In this picture, adding tunneling between the layers turns crossings of site energies into 
avoided crossings of subband energies, and smooths out cusps in the dependences of the 
subband densities on the electric field. 

For double-layer systems, the hopping parameter t is simply related to the subband 
energy separation when the gate voltage is adjusted so that the two layers have equal density: 
2t = E 2 —E\. This simple relationship is very helpful in characterizing experimental systems. 
It is worth remarking that a similar simple relationship exists for triple-layer systems. When 
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Pf = Pb, inversion symmetry guarantees that n\ = n 3 and e\ = e 3 . The Hamiltonian matrix 
can be readily diagonalized for this case with the result: 



E2 = ei, 

£s = + v / ( £i r £i ) 2+2t2 - (3 - 4) 

The subband energies E\ can be determined up to an overall constant by the sublayer 
occupancies N\, which are obtained from the Shubnikov-de Haas (SdH) experiments, using 
Eq. fl2.4| ). Equation ( |3.4| ) can be solved to express €2—61 and t in terms of the two independent 
subband energy differences. We find that 

e 2 -e 1 =E 3 -2E 2 + E 1 , (3.5) 

and that 

(E 3 -E 2 )(E 2 -EJ ^ (3g) 

Since it is clear from Fig. |2| that, for the experimental system, E 3 — E 2 > E 2 — E\ when 
inversion symmetry is established, we see immediately that e 2 — e\ > because of the elec- 
trostatic energy cost of putting electrons in the middle layer. To determine t experimentally, 
it is only necessary to identify the gate voltage V13 at which inversion symmetry is estab- 
lished, and use Eq. ( |3.6|) . V13 may be determined in practice as the point where (Ni — N 2 ) 
is minimized, and electrostatic considerations give V13 ~ V3 + 4ireD-FGN 2 /eo ~ 0.1 V, where 
N 2 ~ 6 x 10 10 cm" 2 is the asymptotic value of N 2 at large Vfg- Unfortunately, the energy 
difference between E\ and E 2 at the symmetric point is close to the limit of resolution of the 
experiment. An estimate of t can be obtained from the minimum difference between N 2 and 
N 3 (for N 3 > 0), which occurs when the densities of layers 1 and 2 are equal. This occurs at 
a front-gate voltage V\ 2 w V s + AireD FG N 3 /e ~ —0.1 V, where N 3 ~ 3 x 10 10 cm -2 is the 
asymptotic value of N 3 at large Vfg- We obtain 

1 . . min(N 2 - N 3 ) , , 
t w -min(E 3 - E 2 ) = — — ~ 0.45 meV (3.7) 

2 2z/q 

for (N 2 - N 3 ) ss 2.5 x 10 10 cm" 2 . 

Subband densities can be extracted experimentally from weak-field SdH oscillation ex- 
periments. The parameters of the model (the offset voltages, the distances to the gates, the 
size-quantization energy e&, and the hopping parameter t) may be determined by fitting to 
the gate-voltage dependence of the measured subband densities. Experimental results are 
compared with a model fit in Fig. |2|. The model calculations were performed with layer 
separation d — 18.4 nm, the mi dwell to midwell distance in the experimental system. The 
offset voltages V^q are consistent with the conditional that the electron density in the layers 
is left-right (LR) symmetric for Vfg — 0.03 V and Vbg — 0, when the total electron density 
is 14.8 x 10 10 cm -2 . Close agreement with experiment is obtained by choosing the tunneling 
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and on-site energies to be within about 25% of t — 0.4 meV and 5% of e& = —4.6 meV. 
The layer densities in Fig. ||] are in close agreement with the SdH data and with the layer 
densities obtained in Ref. [T7| from 3D LDF calculations. This level of agreement provides us 
with the necessary confidence in our model, and determines with some assurance the model 
parameters for the system of immediate interest. We note that this calculation (for B± = 0) 
explicitly neglects nonlocal interlayer exchange and correlation. This approximation appears 
to be well justified at zero magnetic field,@ but we will see that in the strong magnetic-field 
limit, nonlocal interlayer exchange is important. 



IV. QUANTUM-HALL GROUND STATE 

We now turn our attention to the strong magnetic- field limit. Reference finds strong 
QHEs at z/t = 1 and 2; explaining the physics of the underlying incompressible states at 
these filling factors is the main objective of this paper. The fact that a QHE occurs at these 
integer filling factors is, at first sight, surprising. To understand why, it is useful to imagine 
repeating the calculations outlined in the previous section for z/t = 1. The crucial difference 
between the B± = and strong-field situations in such an independent-particle description 
is that, because the kinetic energy is quenched, the density of states consists of a 5 function 
at the subband energies in the strong-field case, whereas it is a constant above the subband 
energy in the B^ = case. As a result, the distribution of density between the three layers 
at v — 1, which we have emphasized is dictated largely by electrostatics, can be expressed 
simply in terms of the subband wave function: 

nj = (27re 2 )- 1 \zf ) \ 2 . (4.1) 

In order to have the charge distributed relatively equally among all the layers, the site 
energies cannot differ by more than ~ t. The difference between subband energies, which 
would give the QH activation gap in such a theory, would then also be ~ t. A gap of this 
size might be reduced or possibly eliminated by disorder in the samples. In the experiments, 
however, the observed gaps can be much larger than t. It seems clear that the explanation for 
these QHEs must lie in the interaction physics of the triple-layer system at strong fields, and 
that the gap would exist even if t were zero. Further evidence is found in the experimental 
observation that as the ratio of side-layer to middle-layer electron density is increased, the 
z/p = 1 state collapses, but the z/p = 2 state state becomes stronger. We shall show that 
these observations are consistent with the behavior expected if spontaneous interlayer phase 
coherence occurs in these triple-layer systems. In the following, we focus on the case z/t = 1; 
the z/p = 2 case is simply related by a particle-hole transformation!^ We shall also assume 
that the system has three layers, although most aspects generalize in an obvious way to 
systems with more than three layers. 

The single Slater-determinant states considered in this section have the form 

( 3 \ 

l*> = n(E¥^j|0>, (4.2) 

where Zj is the component in layer j of a normalized three-subband wave function, and c] x is 
the second-quantization operator that creates an electron in layer j, in the lowest Landau- 
level Landau-gauge state with guiding-center coordinate X. This wave function is a full 
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Landau- level state for the subband with state vector Z^ = (z*, z%, z%). In our HFA, we allow 
Z to be varied to minimize the energy. This variation generally results in a broken-symmetry 
ground state, since z*Zj ^ even when t = 0; i.e., the HF ground state has spontaneous 
interlayer phase coherence. For double-layer systems, the broken symmetry is robust under 
appropriate circumstances, and exists in the exact quantum-mechanical ground state. We 
expect spontaneous interlayer phase coherence to be similarly robust for triple-layer systems. 
It is convenient to define the density matrix 

p jk (X) = (*\c] x Ck X \V) = z*z k , (4.3) 

where |\&) is the coherent ground-state wave function, Eq. (|4.2| ). Note that the filling factor 
of layer j is Uj = pjj. In terms of the density matrix, the HF total energy is 

Erf = J2)-J2 tjkPjk(X) + e jPjj (X) - ^ u a D aj (Y) Pjj (X) 

jX L k aY 

+ W [D jk (X - Y)p kk (Y) Pjj (X) - E jk (X - Y)p kj {Y)p jk {X)]V (4.4) 

Z kY J 

where tj k denotes the tunneling energy between layers j and k, ej is the on-site size- 
quantization energy of layer j, and a indexes neutralizing planes of charges with areal charge 
density ep a = ev a /2n£ 2 produced by remote ionized donors or gates. The unit of length is 



the magnetic length, £ = yhc/eB. As an aside, we note that in the long^wavelength limit, 
Erf has the form of a CP/v-i modeli when expressed in terms of the Zj@ 

In the lowest Landau level, the Coulomb interaction between electrons in layers j and 
k enters through the direct term Dj k (X — Y) and the exchange term Ej k (X — Y). These 
quantities are conveniently expressed in terms of projected Fourier transforms: e.g., 

d 2 q 



D jh (X) = J ±jLD ih {q)e&* J d\(X\v)e^(r\X), (4.5) 
where 



D jk (q) = ^ e -tfU-*l e -4«V (4.6) 



and 



EM = Y J d 2 P D jk (p)e^^ 2 . (4.7) 

The layers are located at z — dj, where j is in general a real number, unless all the layers are 
separated by integer multiples of d. We have neglected the finite thickness of the electron 
wave functions in the z direction, although this effect could be included, if desired. 
When z/ T = 1, it follows from Eq. (pL~^ ) that 



p jk {X)p kj {X) = Pjj {X)p kk {X). (4.8) 

This relation, which follows from the phase- coherent nature of the assumed ground state, is 
of great practical importance to us, because unlike in the B± = case, it allows us to easily 
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express the exchange contribution to the HF ground-state energy E^ F in terms of the layer 
occupancies: 



F,° 

-^hf _ \p 

iv 4> jk L 



-tjk^VjVk + 5jk (tj — ^2 VaD a jj Uj + ~(Dj k — E jk )lJjU k 



(4.9) 



where 

D jk = - 
with v c = e 2 /e £, and 



D (q = 0) - D b -_ fc |(g = 0)] /2nf = -v c ~\j - k\ (4.10) 



= = „ e rdxe-^e-^W*. (4.11) 



v 2vr£ 2 
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Equation (|4.8|) also allows us to write 

p jfc (X) = ^zT-^e^W. (4.12) 

In the ground state, pj k is independent of X, so that we may take e^ (x) = 1. Since 
= 1, the HF ground state is in general determined by two parameters. We will 
concentrate below on the situation where the triple-layer system is symmetric, i.e., z/p = vb- 
Then the system has inversion symmetry around the middle well, v\ = u^, and the ground 
state is completely fixed by z/ 2 . 

The layer occupancies can be calculated by minimizing Eq. ( |4.9|) ; equivalently, the fol- 
lowing simple argument may be used in the absence of tunneling. For sufficiently large, 
all electrons reside in the side layers, closest to the gates. Suppose that both side layers (1 
and 3) are occupied, and imagine moving an electron from layer 1 to the middle layer (2). 
The energy gained by moving to the middle has three components: the on-site energy e^, 
the Hartree energy D^ui + z/ 3 ) — (z/ F + z/ B )] due to the direct Coulomb interaction with the 
side layers and gates, and the exchange energy —E\{v\ + z/3) of the middle electron with 
the side layers. The energy lost by moving to the middle also has three contributions: the 
Hartree energy /^(^ — z?b) of a layer-1 electron with the electrons in layer 3 and the gates, 
the on-site exchange energy —EqV\ in layer 1, and the exchange energy —E^vz of an electron 
in layer 1 with those in layer 3. Specializing to the case of = 1, and using D 2 = 2Dx, the 
energy cost to move an electron to the middle layer is seen to be 

AE = e b + D 1 [{u 1 - z/ 3 ) - (Pp - z/ B )] + (E - E x )v x - {E x - E 2 )u 2 . (4.13) 

The middle layer will be occupied when AE < 0. For the case of equal side-layer densities, 
this happens when e b < e™ ax , where 

e^ = - l -(E,-2E 1 +E 2 ). (4.14) 

A similar argument can be used to find the value e™ m of the on-site energy below which all 
electrons reside in the middle layer: 
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e fn = _(_D 1 -E + E 1 ). (4.15) 



Equations ( [4.14 ) and ( |4. 15| ) show that near- neighbor interlayer exchange (Ei) plays an 



essential role in increasing the size of the interval where all three layers are occupied and 
spontaneous triple-layer phase coherence occurs. The local-density approximation for ex- 
change, commonly used in electronic structure calculations, fails qualitatively for triple-layer 
(and double-layer) systems in the QH regime, because it does not include the effects of in- 
terlayer exchange, included here through Ej for j > 0. The on-site exchange (Eq) favors 
maximizing the charge of individual layers, and E% favors next-nearest-neighbor occupancy. 
For e™ in < ei, < e™ ax , the middle- well occupancy decreases linearly with e&, so that 



(4.16) 



When z/ 2 = 1, all the electrons are in the middle well, and we have, in effect, a single-layer 
system. When i^ 2 = 0, all the electrons are shared between the outside layers, and we have 
a double-layer system. For < z/ 2 < 1, charge exists in all three layers and the HFA ground 
state has triple-layer coherence. 

We may use the HF ground-state calculation of v% in Eq. (|4.16|) to define a phase diagram 
in the space of model parameters in the absence of interlayer tunneling. For symmetric triple- 
layer systems, the state of the system is determined by the middle-well size-quantization 
energy e& in units of the Coulomb energy v c ~ 10 meV, and by the interlayer spacing d 
in units of I. We consider the limit t/v c <C 1, which according to our analysis of the SdH 



data is satisfied by the device studied in Ref. |18|, and set t = 0. The region of stability of 
the triple-layer coherent state at t = is bounded by the dotted and dashed (upper and 
lower) lines of Fig. 0(a). These two lines are defined by the equations z/ 2 = 1 and z/ 2 = 0. 
At fixed d, the system transforms with increasing |e&|, first from a double-layer system to 
a triple-layer system, and finally to a single-layer system with all the charge in the middle 
layer. At fixed e&, these transformations occur in the opposite order with increasing d, as 
electrostatic considerations become more dominant. Figure ^|(b) shows the same plot for 
z/t = 2, obtained from a calculation nearly identical to that of the z/p = 1 case. If, in a given 
sample, front and back gate voltages are adjusted simultaneously to maintain LR symmetry 
as the total density changes, the system will follow a line in this phase diagram that moves 
downward and toward the right with increasing magnetic field, as illustrated in Fig. ^. By 
this procedure, the triple-layer regime should be acessible to experimental study in typical 
TLESs. The HF ground-state calculation by itself suggests that the entire region between 
the dotted and dashed lines in Fig. |^ might support a triple-layer coherent state. However, 
as we show in the following section, a stability analysis using the TDHF equations shows 
that for sufficiently large interlayer spacing, the phase- coherent state cannot be the ground 
state. 

In closing this section, we remark that the distribution of electrons between the three 
layers can be markedly different for the same gate voltages and sample parameters in the 



QH regime, as compared to the B± = case. Reference [18] found, using their 3D LDF 
calculations, that for their triple-layer sample at symmetry, with a total density of 14.8 x 
10 10 cm -2 , the ratio of the density of electrons in the middle layer to the total density was 
z/ 2 /z/t ~ 0.22 when B± = 0. The 2D tight-binding model of the previous section gives the 
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same result. For the same model parameters, we find that in the is? = 1 phase-coherent 
triple-layer QH (3LQH) state, z/ 2 /z/t ~ 0.27, while for z/p = 2 we obtain z^/z^t ~ 0.09. The 
low value of z/ 2 /^t for z/p = 2 is due to the fact that for inversion-symmetric triple-layer 
systems, the A = 2 subband wave function has no weight in the middle layer. Thus B± can 
have a large effect on the ratio of layer densities. 



V. TIME-DEPENDENT HARTREE-FOCK COLLECTIVE MODES 

In the TDHF description of the collective behavior, the variational wave function for the 
ground and excited states of the system both have the HF form, Eq. (|4.2|) , but the density 
matrix elements pj k are allowed to have spatial and time dependence. In particular, the 
low-lying collective modes arise from slow long-wavelength variations of the phase differ- 
ences (pjk(X). In practice, there are several ways to implement the TDHF approximation: 
diagrammatically, by path-integral methods, or by various equation of motion (EOM) ap- 
proaches. TDHF collective modes for a triple-layer system were calculated previously by 
Fertig^l using a diagrammatic approach, for the unphysical case of equidistant, individually 
charge-neutral layers with periodic boundary conditions, with the tunneling and Coulomb 
energies between end layers taken to be the same as between neighboring layers. Here we 
briefly describe a calculation based on the EOM of the density matrix. For the sake of 
generality, we use the language of an iV-layer system, although we will apply our results to 
the case N = 3. 

Following Ref. [H], we define the projected density-matrix operator 



PM=eW£$x [J d 2 r(X\r)e-^(r\Y) 



XY 

and the density-density response function 



CkY, (5.1) 



XjMm(<l,t) = *(Tpjk(q,t)pjk(-q,0)), (5.2) 

where T denotes the time-ordering operator. The physical quantities of interest to us are a 
function of q — |q|, due to the isotropy of the system. The HFA to \ is obtained from the 
EOM 

ihd t p jk (q) = [p jk (q),Hn F ], (5.3) 
where is the HF Hamiltonian. The commutator is evaluated by use of the identity 

[p jk (pl Pzm(q)] = hiPjmiP + q)e^ px ^ 2 - 6 mj p lk {q+ p)e^ qx ^ 2 . (5.4) 

In the TDHFA, x(g, ui) is obtained from the HF x by including the effects of the 
Coulomb interaction between a particle in layer j and a hole in layer k with total in-plane 
momentum q, within a generalized random-phase approximation. This results in a Dyson 
equation for \ in terms of the direct and exchange interaction between a particle and hole. 
In matrix notation, 

X = X (0) + X {0) \w X , (5.5) 
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where W = H — X is the sum of the Hartree and Fock contributions to the particle-hole 
interaction, given by 

Hjkim{q) = S jk 5i m D ^J , (5.6) 

and 

Xjkim(q) = 5jm,8ki • (5.7) 

Since the layer indices each have N possible values, Eq. ( |5 . 5|) may be solved numerically by 
inverting an N 2 x N 2 matrix. 

The collective modes are obtained from the poles of x{q, u ) that have nonzero residues. 
The results of a sample calculation for (t = 0,—€b/v c = 0.75, d/t = 1.6), which has 
v-i = 0.838, are shown in Fig. |j. (The collective mode for = 2 can be obtained from 
that for z/x = 1 by particle-hole conjugation, using u 2 — > 1 — z^-) As seen from Fig. 01, there 
are two collective modes, which correspond to linear superpositions of variations in 12 (A) 
and 023 P0- As the spacing d between the layers increases, one of the collective modes 
softens and eventually becomes unstable at a wavevector q ~ l~ x . This signals the onset of 
a charge- density wave instability in the HFA, on which we comment further below. The lack 
of an anticrossing repulsion at the point (away from q — 0) where the two collective-mode 
frequencies are equal is a special feature of the assumed LR symmetry of the charge distri- 
bution in the TLES. When this symmetry is present, the collective modes are excitations 
from a subband state (A = 1) that has even parity to subband states that have even (A = 2) 
or odd (A = 3) parity. The opposite parity between the two final states provides a selection 
rule that prevents the two collective modes from mixing. 

So far, we have focused on the case where the tunneling tjk is negligible compared to the 
Coulomb interaction energy scale v c , and can be set to zero. This results in gapless Goldstone 
modes, as shown, for example, in Fig. |. (When 1 < i/ T < N — 1, gapped collective modes 
exist, even in the absence of tunneling.) The tunneling energy can be varied by changing 
the thickness of the barrier between the quantum wells. The [^(l)]^ -1 invariance associated 
with the freedom to choose the N — 1 relative phases <pj,j+i in each layer is broken once the 
electronic states in different layers are coupled by tunneling. The loss of this invariance gives 
a gap of order 2t (for nearest-neighbor tjj. = t) in the collective- mode spectrum as q — > 0. 

In our TDHF collective-mode calculations, we find that for fixed e&, the collective-mode 
spectrum in the triple-layer regime softens with increasing d, and that, except for sufficiently 
small values of d/£, an instability occurs before the z/ 2 = line in Fig. |3](a) is reached. 
The first collective mode that goes soft corresponds to an excitation from the filled A = 1 
subband to the empty A = 2 subband, which has all the charge on the outside layers; this is 
the favored arrangement for the charge at large d. This instability appears in the TDHFA as 
an imaginary-valued collective-mode frequency. Based on extensive calculations performed 
previously for the case of double-layer systems,^ the broken translational symmetry HF 
ground states that are found for larger values of d will quickly lose their interlayer phase 
coherence, and the charge gap (incompressibility) necessary for the integer QHE will rapidly 
go to zero. In our view, the broken translational symmetry of the state on the large-c? side of 
this instability is likely to be an artifact of the HFA, which can enhance intralayer correlations 
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only by breaking translational symmetry. When quantum fluctuations are included, the 
broken translational symmetry of this state is likely to be lost, but in our view the loss of 
interlayer phase coherence and the vanishing of the QH charge gap will remain. We therefore 
use the location of the TDHF instabilities as an estimate of the layer separation at which 
spontaneous coherence and the integer QHE are lost. The results are shown in Fig. |^(a) 
for z/p = 1. The solid line shows the border of the TDHF instability. For z/p — 1, the 
region between the upper dotted line (^=1) where the middle layer is fully occupied, and 
the solid line, where the TDHF instability occurs, is the resulting estimate of the region in 
parameter space where a 3LQH coherent state occurs. This shows that as djl increases, the 
z/p = 1 phase- coherent 3LQH state is more likely to be stable when the middle layer has 
more electrons than the side layers.0 Figure §(b) shows the phase diagram for z/p — 2. In 
this case, the 3LQH state at larger d/£ is likely to be more stable when the middle layer 
has fewer electrons than the side layers. These behaviors are seen in the experiments of 
Ref. [H| The dot-dashed lines in Fig. [5] represent values of (d/£, — £&/u c ) for a hypothetical 
sample with d = 18 nm, as the total density is varied from 1 to 2 x 10 11 cm -2 , where we 
have used 6b = —8 meV for z/p — 1 and e& = —6 meV for z/p = 2. As the total density 
and B± are increased while keeping the total filling factor z/p constant, the dot-dashed lines 
in Fig. |3| show that the ratio of middle-layer density to side-layer density decreases. In 
all cases, the presence of interlayer tunneling will enlarge the region where a triple-layer 
QHE is expected to occur. We emphasize that the TDHF instability with increasing d is 
driven by the increasing relative importance of intralayer correlations compared to interlayer 
correlations. At large enough layer separations, the interlayer coherence that provides for 
good interlayer correlations will be lost, along with the charge gap necessary for the QHE. 
This scenario for the disappearance of the QHE with increasing layer separation does not 
require that the TLES make a transition to a bilayer system, as was hypothesized in Ref. ITS. 



VI. CONCLUDING REMARKS 

We have shown that TLESs, such as those fabricated by Shayegan and co-workers, can 
be described using a simple tight-binding model, which allows them to be characterized in 
terms of a small number of parameters. The tight-binding model is able to quantitatively 
account for weak-field SdH experimental results for the dependence of the three subband 
energies on the gate voltage. Using this model, we have estimated the dependence of the 
ground state of the system on model parameters for filling factors z/p = 1 and 2. We propose 
that in triple-layer systems, as in double-layer systems, the QHE can occur at integer total 
filling factors, even when the charge is distributed among all three layers and there is no 
hopping between the layers. (Tunneling would be required for noninteracting electrons to 
produce a gap at integer filling factors.) These QHEs occur because of the formation of 
broken-symmetry ground states with spontaneous interlayer coherence. 

The distribution of charge between the three layers of the system is typically determined 
predominantly by electrostatic considerations, and states can occur with the electrons dis- 
tributed among one, two, or three layers. Our HF theory of the ground state shows that 
for physically accessible regions of the middle-well on-site energy and layer separation, a 
triple-layer phase coherent state is possible. The stability of this state was estimated using 
the TDHFA, and a phase diagram was constructed, delimiting the regions in the parameter 
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space of the TLES for which triple-layer coherence is likely to be found. The QH states 
studied by Shayegan and co-workers at z/p = 1 and 2 show the behavior expected for phase- 
coherent states. In particular, they exhibit the QHE, even for small tunneling energies and 
unequal layer densities. Additionally, the z/ T = 1 QHE is suppressed by increasing occu- 
pancy of the side wells at the expense of the middle well, while for z/p = 2, the opposite 
is true, in agreement with our findings. Future detailed comparisons between our phase 
diagram and experiment will be facilitated by the possibility of following lines in the phase 
diagram for a single sample, by adjusting gate voltages so that inversion-symmetric Ut — 1 
and 2 states occur for a range of magnetic fields. In any such quantitative comparison, it 
will be necessary to approximately account for the finite widths of the individual quantum 
wells, which we have not done here. 

As in the case of the double-layer coherent state, a useful footprint of interlayer coherence 
is unusual sensitivity to small tilts of the magnetic field away from the normal to the layers.0 
The most convincing evidence for triple-layer coherence would be the observation of a strong 
suppression in the activation energy of the 3LQH states due to the application of a moderate 
parallel magnetic field. The state of the triple-layer system in a tilted magnetic field is 
most informatively described using a field-theoretical approach. We have calculated the 
magnitude of the parallel field required to suppress the activation energy as a function of 
the layer density, for both triple-layer and unbalanced double-layer systems. The effect 
of a parallel field on triple-layer states, charged order-parameter textures of triple-layer 
states, and other related properties of triple-layer phase-coherent states will be discussed in 
a forthcoming paper.il 
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FIG. 1. Schematic illustration of the triple-layer electron system of Ref. 18. The solid lines 
represent the energy of the confining barriers, and the long-dashed lines are the energies of the 
lowest-energy quantum state for a given well. The dotted curves represent electron densities, which 
are peaked at midwell. Our model idealizes the electron density to be midwell 5 functions. 
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FIG. 3. (a) Phase diagram for the symmetric i>t = 1 3LQH coherent state, (b) Phase diagram 
for the symmetric is? = 2 3LQH coherent state. The mean-field regions of stability lie between 
the dotted and solid lines. All calculations are for t = 0, for the case of equal left and right 
layer densities. The dot-dashed lines represent values of (d/£,—eb/v c ) obtained for a sample with 
d = 18 nm, for total density in the range of 1 to 2 x 10 11 cm" 2 , using et = —8 meV for zv T = 1, 
and efe = —6 meV for = 2. 
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FIG. 4. Collective modes of the v>t = 1 3LQH coherent state for (t = 0,—eb/v c = 0.75, 
d/i = 1.6), which has = 0.838. The lack of an avoided crossing at ql ~ 0.4 is due to the 
assumed left-right inversion symmetry of the system. 
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